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RANDOM  FUNCTIONS  FROM  A 
POISSON  PROCESS 

ROBERT  FORTET 

UNIVERSITE  DE  CAEN 

# 

1.  Definition  of  random  functions  from  a  Poisson  process 

Let  R(t,  r )  be  an  ordinary  (nonrandom)  function  of  two  variables  t  and  t, 
which,  for  example,  might  represent  time.  Let  us  consider  instants  distributed  at 
random  on  the  time  axis  according  to  the  classical  Poisson  process  of  density 
m  (m  >  0).  Let  iVi(0  be  the  number  of  these  instants  belonging  to  the  interval 
(0,  t)  if  t  >  0,  to  the  interval  (/,  0)  if  /  <  0.  We  shall  write  N(t)  —  N\(t)  if  t  >  0; 
N(t )  =  —Ni(t)  if  t  <  0,  so  that  N(t')  —  N(t)  represents  the  number  of  instants 
belonging  to  ( t ,  t')  whatever  may  be  t  and  i'  (t  <  t'). 

The  following  well  known  results  should  be  recalled: 

(a)  N (t)  is  a  random  function  with  independent  increments,  taking  integer  values, 
and  is  almost  certainly  nondecreasing. 

(b)  The  probability  that  the  number  of  instants  belonging  to  any  finite  interval 
of  time  be  infinite  is  zero. 

(c)  The  distribution  in  the  time  of  these  instants  is  stationary. 

(d)  m  is  the  expectation  of  the  number  of  instants  belonging  to  any  interval  of 
amplitude  1. 

(e)  The  probability  that  the  number  of  instants  belonging  to  any  interval  (/,  i'), 
t  <  t',  be  equal  to  n  is  equal  to 

a  (jn8)^ 
n\  ’ 

where  8  =  t'  —  t. 

Calling  n,  T2,  ,  Tj,  .  .  .  the  instants  of  an  interval  (to,  t),  (t  >  r0),  different 

applications,  for  example,  the  phenomenon  of  noise  in  electronics,  lead  to  a  con¬ 
sideration  of  random  functions  of  t  defined  in  the  following  way: 

(1)  X(t) 

which  may  be  written  as 

(2)  X(t) 

More  generally,  we  shall  consider 

(3)  X(t)  =  f+a>R(t,  r)dN(r), 

and  we  shall  denote  random  functions  of  this  kind,  that  is  random  functions  from 
a  Poisson  process,  by  P.r.f.  The  study  of  these  functions  has  been  developed  by 


=  ^R(t,  T,'), 
i 

=  Jr (t,  T)dN(r) 
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several  authors,  in  particular,  Blanc  Lapierre  [1],  Kac  and  Hurwitz  [3]  and  Maru- 
yama  [2].  But  it  is  possible  to  extend,  to  simplify  and  to  complete  their  results;  this 
is  the  purpose  of  our  paper. 

The  first  point  is  to  give  a  precise  definition  of  integrals  (2)  or  (3).  They  exist 
only  if  R(t,  r)  satisfies  some  conditions  which  depend  on  the  sense  (integral  with 
probability  one,  integrals  in  quadratic  mean,  .  .  .)  attached  to  the  integrals  (2),  (3). 
Let  R(t )  be  a  function  of  r.  We  first  consider  the  integrals 

(4)  Xa0=  f*R(T)dN(T), 

•'a 

where  a  and  (a  <  0)  are  any  finite,  given  constants.  Denote  by  assumption  Hi 
the  assumption  that  R(t)  is  a  definite  function  which  is  .^-measurable  on  (—  00 , 
+  00 )  and  finite,  except  perhaps  on  a  set  of  measure  zero. 

1.1  Integrals  with  probability  one.  We  can  give  to  (4),  considered  as  an  integral 
with  probability  one,  two  different  meanings. 

(a)  Except  in  cases  whose  probability  measure  is  zero  any  trial  e  assigns  to  N  a 
finite  number  n  of  jumps  (n,  .  .  .  ,  r„)  between  a  and  /3.  It  is  known  that  if 
2V(/3)  —  N(a )  =  n,  the  tj  are  independent  and  follow  the  uniform  law  on  (a,  /3). 

We  can  consider  the  sum  R(ry),  and  it  follows  from  assumption  Hi  that  with 

i 

probability  one  this  sum  exists  and  is  finite.  Then  the  formula 

(5)  X+-  tj)  for  e 

J 

defines  a  random  variable  that  we  may  consider  to  be,  by  definition,  Xap. 

(b)  We  can  define  X^  as  the  stochastic  limit,  if  it  exists,  of  a  Stieltjes-Riemann 

k 

sum  of  the  following  type:  ^  j?  ( £,)  A SN,  obtained  by  dividing  (a,  /3)  in  partial 

*— 1 

intervals  by  points  t0  =  a,  t8+i,  .  .  .  ,  <  4  =  /3,  putting  A „N  = 

N(tg )  —  N(tg- 1),  £s  being  chosen  arbitrarily  (not  at  random)  on  (/„_ i,  is).1  The 
limit  is  obtained  when  sup  (/,  —  4_i)  — >  0.  It  is  necessary  that  this  limit  exist  and 
be  independent  of  the  division  points  ( t„ ),  and  of  the  choice  of  the  £*’s.  Definition 

(а)  is  more  suitable  to  the  physical  origin  of  the  considered  integral.  Definition  (b) 
is  more  suitable  to  the  usual  definition  of  stochastic  Riemann,  or  Stieltjes-Riemann 
integrals,  but  it  requires  an  additional  condition  on  ./?(£*),  in  order  that  the  inte¬ 
gral  exist  with  probability  one ;  indeed  for  a  trial  e  such  that  N(f3)  —  AT(a )  =  n, 
if  sup  (/*  —  4_i)  — » 0,  there  will  be  a  time  when  (5)  will  be  confused  with 

(6)  2>(w, 

j 

&  denoting  a  number  which  tends  to  ry,  but  in  an  arbitrary  manner.  In  order  that 

(б)  has,  with  probability  one,  a  limit  independent  of  the  choice  of  the  £s’s,  namely 

R(tj),  it  is  necessary  and  sufficient  that  the  set  of  points  of  (a,  0)  where  R(t) 
i 

is  discontinuous  be  of  measure  zero.  Let  H 2  denote  the  assumption  that  the  points 
1  We  may  consider  choosing  the  £’s  at  random. 
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of  (—  <» ,  +  oo )  at  which  R(t)  is  not  continuous  form  a  set  of  measure  zero  (Hz  im¬ 
plies  Hi). 

It  is  seen  that  whatever  may  be  a  and  /S,  the  integral  with  probability  one 

(4)  exists  in  the  sense  (a)  if  Hi  holds,  and  in  the  sense  (b)  if  Hz  holds.  Now  let  Hz 

/+00 

|  R(t)  |  dr  <  -f-  oo  (this  integral  is  a  Lebesgue 

■  OO 

integral).  Let  A  be  a  random  variable  following  the  uniform  law  on  (a/3),  and 
consider  the  random  variable  J?(^4).  If  Hz  holds,  the  expectation  of  R(A)  exists 

J/'P 

'  R(t)<It.  It  is  seen  that  when  N(f. 3)  —  N( a)  =  n, 

a 

the  conditional  expectation  of  Xap  exists  and  is  nn/i 3  —  a;  then,  the  a  priori  ex¬ 
pectation  of  Xap  exists,  and  is  m/x.  This  is  true  with  definition  (a)  as  well  as  with 
definition  (b). 

Now  let  X  be  the  integral  with  probability  one, 

(7)  X  =  f+C°R  ( r)  dN  (  t)  . 

j  —  00 


We  will  define  it  as  the  limit  with  probability  one,  if  it  exists,  of  X ^  as  a  tends 
to  —  oo  and  /8  tends  to  +  oo  independently.  It  follows  from  a  classical  theorem 
that  X  exists  in  this  sense  if  Hi  and  H 3  (or  Hz  and  Hs)  hold.  It  is  seen  that  the 
expectation  of  X  exists  and  is  the  limit  of  the  expectation  of  Xap,  so  it  is  equal  to 

/*  +00 

ml  R(t)  dr  . 


1.2.  Probability  law  of  X^  and  of  X.  In  the  case  in  which  assumption  Hi  (or  Hz) 
holds,  let  a  be  a  given  number  and  consider  X^  as  a  function  of  the  variable  /8  (it 
is  obvious  that  we  can  exchange  the  role  of  a  and  j8).  Obviously  X^  is  a  continuous 
in  probability  random  function  of  /8  with  independent  increments;  indeed,  if  we 
write  AN  —  N(fi  +  A/3)  —  iV(/3),  AX  =  Xa,i 3+49  —  X^,  we  get  AX  =  0  as 
AN  =  0;  but  the  probability  of  AN  =  0  tends  to  1  as  A/8  0.  We  then  have  as  a 

result,  that  the  probability  law  of  Xap  is  an  infinitely  divisible  law  and  then  the  prob¬ 
ability  law  of  its  limit  X  is  also  an  infinitely  divisible  law. 

Let  6(v)  be  the  characteristic  function  of  R(A).  The  conditional  characteristic 
function  of  Xo/j  as  N($)  —  N( a)  =  n  is  obviously  [0(d)]".  Then  the  a  priori  char¬ 
acteristic  function  of  Xap  is 

^  g-m(P--a)  m.  {fL. — OlL  [  0  (  v)  ] n  =  exp  { m  (/3  —  a)  [  6  (  v)  —  1  ]  } . 

n 


The  logarithm  pap(v)  of  the  characteristic  function  of  X<^,  is 

(8)  ypaai  v)  =  m  (/8  —  a)  [6  (  »)  —  1]  =  m  f  [  eivRW  —  1]  dr  . 


From  this  there  results,  and  it  is  easy  to  verify  directly,  that  if  not  only  Hi  (or  Hz) 
holds,  but  also  Hz  holds,  the  integral 


eivR(r)  _  1]  dT 
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converges  as  a  and  /S  tend  independently  to  —  and  +  °°  respectively.  Hence, 
the  logarithm  \p(v)  of  the  characteristic  function  of  X  is  equal  to 


(9) 


\p  (  v)  =  m  J  [  eivR(T~)  —  1  ]  d  t  . 


These  results  will  be  useful.  We  see  that  they  permit  the  computation  of  the  mo¬ 
ments  of  Xafi  and  of  X,  when  they  exist.  In  particular  the  “cumulants”  of  X,  when 
they  exist,  are  easily  computed  from  the  integrals 


d  t  . 


1.3.  First  generalization.  Let  Ri(t),  i^2(r),  .  .  .  ,  Rk(r)  be  k  functions  such  that 
Hi  (or  Hi)  and  Hz  holds.  The  k  random  variables 

/+OO 

Rj(r)  dN  (  r)  (integral  with  probability  one) 

.00 


are  not  independent.  Let  if/(v  1,  ,  vk)  be  the  logarithm  of  the  characteristic 

function  of  the  ^-dimensional  distribution  of  the  ^-dimensional  random  variable 
Z  —  (Xx,  .  .  .  ,  Xfc).  Let  us  set 

Y  =  V\X\  +  ...-+-  vkXk  . 


We  have  ip(vi,  .  .  .  ,  vk)  =  log  {E[elY}). 

/•  -bco  /» -boo 

F=  /  (»,J?,(t)  +  .  .  .  +  vtRk(.r)]iN(r)  =  /  R(r)dN(r), 

J CD  */  —  00 

where  R(t)  =  ViRi(t)  +  .  .  .  +  VkRk(r). 

For  any  system  (*>i,  .  .  .  ,  vk)  with  v/s  finite,  R(t)  satisfies  Hi  (or  Hi)  and  Hz’, 
hence,  from  (9) 

(9')  t(vi ,  ...  ,vk)  =m  f+C°  [e^  -  1]  dr. 

This  expression  can  be  used  to  compute  the  moments  of  Z. 

1.4.  Second  generalization.  In  the  integrals  (4)  or  (7)  we  ascribe,  to  a  jump  of 
N(t)  at  the  abscissa  t,  a  given  extent  R(t).  We  may  imagine  that  this  extent  is 
chosen  at  random;  in  other  words  we  may  attach  to  each  value  of  r  a  random 
variable  Y (r)  and  consider  the  integrals 

X^  =  £>YI,T)dN(r),  X'  =  f_+jY(T)dN(r). 

Xi#  may  be  defined  (definition  a)  as  equal  to  F (r})  for  a  trial  e,  in  which  the 

i 

jumps  of  N{t)  between  a  and  0  occur  at  the  points  n,  .  .  .  ,  17,  .  .  .  ,  and  X'  may 
be  defined  as  the  stochastic  limit  of  X'a0  as  a->  -  00,  (3— >  +  « .  If  the  F (r)’s  are 
mutually  independent,  the  law  of  X^  will  be  again  an  infinitely  divisible  law,  and 
also  the  law  of  X';  besides  it  may  be  sufficient  that  F (r)  be  defined  except  on  a  set 
of  values  of  measure  zero.  Let  us  make  the  following  assumptions:  except  perhaps 
on  a  set  of  values  of  measure  zero,  F(r)  is  defined ;  £[  |  F(r)  |  ]  and  E[  Y (r)]  exist  and 
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are  ^-measurable  functions  of  r;  X'  exists  as  the  limit  almost  certainly  of  and 

/+00 

E[Y(r)]dT. 

00 

Let  us  call  6(t;v)  the  characteristic  function  of  F(r).  Assuming  that,  for  any  v, 
6(t;  v )  is  an  ^-measurable  function  of  t,  we  obtain  for  the  logarithm  yf/'(v)  of  the 
characteristic  function  of  X' 

(10)  t'(v)  —  m  f  — 1  ]dr. 

J  —  CO 

As  an  example  let  F(u)  be  a  distribution  function  such  that 

/*  +00 

udF ( u )  = X . 

that  the  existence  of  the  expectation  X 
as  a  — >  +  00  • 

1  -  F  {t  -  r)  , 

F  (/  -  r) 

when  t  <  t,  and  F(r)  =  0  almost  certainly  when  r  ^  t;  for  example  X'  might 
represent  the  number  of  conversations  which  are  being  conducted  at  time  t  in  a 
telephone  center,  if  it  had  no  loss.  Let  F(u)  be  the  distribution  of  the  length  of  the 
conversations.  We  have 

6  (r,  v )  —  eiv  +  F  (t  —  r)  [1  —  eiv]  if  t  <  t 

6  (r,  v)  =  1  if  r  ^  . 

It  is  easy  to  verify  that  the  above  assumptions  hold,  and  that 

yj/'  ( v )  =  m\  [eiv  —  1] , 

so  that  the  number  of  conversations  conducted  follows  a  Poisson  law. 

1.5.  Integrals  in  quadratic  mean.  We  may  also  consider  the  integrals  (4)  and  (7) 
as  integrals  in  quadratic  mean.  We.  shall  call  H\  the  assumption  that  R(t )  and 
R2(t)  are  integrable  in  the  sense  of  Riemann  on  (—  °° ,  +<»).  Suppose  that  Hi 
holds  (this  implies  that  Hi,  H 2  and  Hz  hold),  and  let  us  consider  the  Stieltjes-Rie- 
mann  integral  (5).  To  prove  that  this  integral  converges  in  quadratic  mean,  we 
must  prove  that 

£[j2  A. at}2], 

9 

where  £»,  are  arbitrary  numbers  on  (/,_  1,  t,),  tends  to  zero  as  the  supremum  of 
(/«  —  4-i)  tends  to  zero. 

If  we  write  us  =  i?(£g)  —  !?(£'),  the  considered  quantity  may  be  written  as 

(4  —  4- 1)  +  »»2[  2  u*  (*»  —  *.-1)  ]  —  (4  “  *»-i) 2] 


F( 0)  =  0  ;  / 

"'o 

We  have  f  udF(u)  >  a[l  —  F(a)j,  so 

•'a 

implies  that 

a  [1  —  F  (a)  ]  — >  0 
Suppose  that  F (r)  be  such  that 

Pr  [Y  (r)  =  1]  = 
Pr  [Y  (r)  =  0]  = 
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which,  under  assumption  Hi,  tends  to  zero.  We  see  also  that 

Js[j  ^R(^)AaN\2]-^mfaPRHT)dT-{-m2[f^R(T)dTj. 

Then  this  last  quantity  is  the  second  moment  of  X^.  It  follows  that  if  a  — ■»  —  <» 
and  0  — >  +  00  (independently  each  other),  tends  in  quadratic  mean  to  X,  the 
second  moment  of  which  is 

(11)  m  J'+C°R2  (  t)  d  t  -f-  J^+C°R  (t)  d  . 

Remark.  It  may  be  noted  that  the  continuity  of  R(t)  is  not  necessary  to  obtain 
the  above  and  following  results;  nevertheless,  we  know  that  a  Riemann  integrable 
function  is  continuous  almost  everywhere. 

1.6.  Another  point  of  view.  Set 

N*  (/)  =  N  (t)  -  EN  (t)  =  N  (t)  -  ml . 

In  many  applications  we  have  to  consider  not  Xap  and  X,  but  the  quantities 
(4')  X*B=  fPR(T)dN*(T )  and  X*  =  f+<°R  (  r)  dN*  (  r) . 

The  integration  in  quadratic  mean  leads  us  to  a  study  of 

£[(  2  =  m  2  u\  (/.  -  tg—\) . 

If  we  call  H6,  which  is  much  less  restrictive  than  Hi,  the  assumption  that  R2(r)  is 
integrable  in  the  sense  of  Riemann  on  (—  °° ,  +  °° ),  it  is  easy  to  see  that  a  sufficient 
and  also  necessary  condition  for  the  existence  in  quadratic  mean  of  X*  is  that 

+00 

(110  E[X*]  =0,  E[X**]=m  R2  (r)  dr . 

J — 00 

/+00  f 

R2(r)dT  exists  in  the  sense  of  Lebesgue  (assumption  Hh, 

•  00 

less  restrictive  than  Hf).  X^  and  then  X*p  exist  almost  certainly  [in  the  sense 
of  (a),  p.  2],  and  we  have 

£[xy=°,  Eix»]  =mffRHr)dr. 

Then  X*  exists  as  the  limit  in  quadratic  mean  of  X*fi,  and  (110  has  a  meaning. 

1.7.  Conclusion.  Obviously,  it  is  sufficient  to  apply  the  above  results  to  solve 
the  problem  of  the  definition  and  of  the  existence  of  the  P.r.f.  X(l)  defined  by  (3). 

2  Convergence  to  a  Laplace  law 

Suppose  that  m  —>  +  <» .  What  are  the  conditions  under  which  the  law  of  X  con¬ 
verges  to  a  Laplace  law?  That  is  to  say,  what  are  the  conditions  under  which  there 
exists  a  number  a,  and  a  positive  number  b  (depending  on  m)  such  that 

lim  [p  (  v/  b )  —iav/  b]  —  —  v2/  2  . 

m — >+oo 


(12) 
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p(v)  is  given  by  (9).  Suppose  that  H 4  holds.  Differentiating  under  the  integral 
sign,  we  get 

\f/'  (  v)  =  im  f  R(  t)  elvR^d r ,  (0)  =  im  f  R(t)  dr  , 

*/—co  */— 00 

V  (  v)  =  -m  f+a>RHr)  eivR^dr,  ^"(0)  =  -m  f+<°R*(T)dT. 

J — CO  J _ CO 


We  note  that  the  differentiations  are  justified,  that  p"(v)  is  a  uniformly  continuous 
function  of  v. 

Setting 

Mi  =  f  R(T)dr,  M2=  f  i?2(r)dr, 

we  get 

fft 

fi(v)  =  im/iiv  —  -y  M2^2+  v2u  (  v) , 


where  co(z»)  — »  0  as  v  — >  0.  Excluding  the  noninteresting  case  ju2  =  0,  we  see  that 
(12)  holds  when 

a  —  mu  1  ;  b  =  s/mm  . 


Remark.  It  is  easy  to  extend  this  result  to  the  case  considered  in  the  second  gen¬ 
eralization  above,  but  we  shall  consider  the  case  of  the  first  generalization.  As¬ 
suming  that  H4  holds  for  i?i(r),  .  .  .  ,  Rifr),  the  same  method  shows  that  if 


/+00  S'  -TOO 

Rj(r)  dr  ,  Mj.«=  /  Ry(T)R,(r)  dr  , 

.00  •/  —  00 

. _  PI  n  ■ 

the  ^-dimensional  random  variable  made  up  of  the  random  variables  — 

tends  as  m—>  +  00,  to  a  Laplace  variable;  indeed,  the  logarithm  of  the  charac¬ 
teristic  function  of  its  limit  distribution  is 


1 

2 


)■ 


^2  HjltVjV,  is  a  definite  positive  quadratic  form  of  the  vfs  if  the  i?/r)’s  are  linearly 
j.* 


independent,  because 


f_m  y^Ri(T)  yi]  dr- 


2.1.  Another  point  of  view.  We  may  have  the  convergence  to  a  Laplace  law  under 
different  conditions.  Xap  being  a  random  function  with  independent  increments, 
may  be  considered  as  the  sum  of  independent  random  variables.  Then,  as  a— ■»  —  °°  , 
/3  — ♦  +  00 ,  it  will  be  asymptotically  a  Laplace  variable,  at  least  under  some  condi¬ 
tions  which  may  be  easily  found  from  the  general  theory  of  the  convergence  of  a 
sum  of  independent  random  variables  to  a  Laplace  variable.  For  example,  let 
Ri(t)  .  .  .  Ric(t)  be  k  functions  such  that  Rfir )  and  R)(t)  have  finite  Riemann  in¬ 
tegrals  on  any  finite  interval  (a,  ($),  (j  =  1,  2,  .  .  .  ,  k),  but 

J  R)(t)  dr  =  +  00  . 
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Let 
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=  dr  ,  (  r)i?,  (  t)  r  , 

and  suppose  that  the  limits 


lim 


Mi>» 


llXXi  /— - ~ 

a — >  —  00, S — >  +  ao  V  JU  .  .  U 

~  'j ,  J  a, a 


=  Mi.- 


(j,  s  =  1  ,  .  .  .  ,k) 


exist  (clearly  they  exist  and  take  the  value  1  if  j  =  s). 
Let  be  the  subset  of  (a,  0)  on  which 

l-Rj  (T)  I  >e 

and  suppose  that  whatever  may  be  j  and  e,  we  have 


lim 

a — >  —  00  ,/3 — >  +  00 


=  0 


From  (9')  it  is  easy  to  show  that  the  ^-dimensional  random  variable  composed  of 
the  set  of  random  variables 

f  R)(t)  dN  (r)  —  mu'j 

*'a 

y/1^- 

'jj 

has  a  distribution  which  tends,  as  a  — *  — >  +  °o,toa  yfe-dimensional  Laplace 
distribution  whose  logarithm  has  the  characteristic  function 


A  result  of  the  same  kind  may  be  given  when  one  of  the  limits  a,  0  is  a  constant 
and  the  other  tends  to  infinity. 

Remark.  Above  and  below  we  suppose  that  the  functions  R(t,  r),  R(t),  Rj(t) 
are  real,  but  it  is  not  difficult  to  extend  the  results  to  the  case  where  R(t),  Rj(t), 
R(t ,  r)  are  complex  functions. 

2.2.  Application .  It  is  sufficient  to  apply  the  above  results  to  get  the  main  re¬ 
sult  concerning  the  P.  random  functions  found  by  Kac  and  Hurwitz,  Blanc  La- 
pierre,  and  Maruyama.  For  example,  let  us  consider  the  random  function  X(t) 
defined  by  (4'),  and  suppose  that,  for  any  t,  R(t,  r)  considered  as  a  function  of  r 
verifies  the  assumption  E 5. 

X*(t)  =  f+C°R  (t,  r)dN*(r) 

J  —  00 


defines  a  random  function  of  t.  The  expectation  of  X*(t)  exists  and  is  zero. 
Its  variance  is 

<r2  [**(*)]  =  m  f  °°i?2(/,  r)  d  t  =  m<r2  (t) , 

J  —  00 


where 
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It  follows  that  X*(t)  is  a  random  function  of  second  order;  its  covariance  T(/,  t')  is 
found  from  (9')  by  setting  k  =  2,  Ri(t)  —  Ri(t,  r),  Rz(t)  =  R(t',  r)  in  (9').  Per¬ 
forming  the  substitution,  we  have 


where 


T  (/,  t')  =  m  f  R  ( t ,  t)R  ( t ',  t)  dr  =  mmit,  t') , 

J  —  CO 

/*  +00 

M2  a  t')  =  R(t,  r)R(t',  t)  dr  . 

J  —  00 


In  the  particular  case  (very  important  in  the  theory  of  noise)  of  (2),  where 
R(t,  t )  =  R(t  —  t)  is  a  function  only  of  /  —  r,  with  R(t,  r)  =  0  for  r  <  i0  or 
t  >  l,  and  R(t,  r)  =  R(t  —  r)  for  to  <  t  <  t,  assuming  t'  >  t>  to,  we  get  from  (13) 

M2(U')  =  flR(t—  t)R  (tf  —  t)  d  t  =  °R  (u)R  (/'—/+  u)  du  . 

*'1«  */o 

For  to  =  —  00 , 

M?  (^>  t')  =  f  R(u)R(t'  —  t  +  u)  du  , 

•/—CO 

which  depends  only  on  t'  —  t,  and  in  fact  only  on  |  /'  —  / 1 .  Besides,  under  the  same 

/•  +00 

conditions  <r2(/)  is  independent  of  /.  This  implies  that  X (t)  =  /  R(t  —  r)dN(r) 

is  a  stationary  random  function  of  the  second  order.  This  was  obviou^and  further 
it  is  obvious  that  X*(t )  is,  more  precisely,  a  strictly  stationary  function. 

In  practice  m  is  generally  very  large,  and  we  saw  that,  as  ra— >  -f  ,  the  tem¬ 
poral  law  of  the  process  X*(t),  suitably  normed,  tends  to  the  law  of  a  Laplace 
process,  the  correlation  function  of  which  is 

/*  +00 

(14)  r  (h)  —  /  R  (u)R  («-J-  h)  du  . 

So  for  any  function  R(u)  satisfying  H'h,  the  integral  (14)  is  a  positive  definite  func¬ 
tion  of  h.  It  is  known  also  that  r(h )  is  a  continuous  function  of  h;  in  other  words, 
r(h)  is  a  characteristic  function.  Conversely,  given  any  characteristic  function  r(h), 
we  would  like  to  find  a  function  R(u)  so  that  (14)  holds.  From  a  classical  theo¬ 
rem  [5,  p.  91]  it  is  necessary  that  the  spectral  function  F( co)  of  r(h)  has  a  density 
function /(co)  (be  absolutely  continuous).  In  this  case  R(u)  exists,  and  is  defined  as 
the  Fourier  transform  of 

g( co)  =  V/  (co) . 

R(u )  obtained  in  that  way  belongs  to  L2,  that  is  to  say,  satisfies  H'h.  Then  we  can 
say  that  every  Laplace  stationary  process  with  absolutely  continuous  spectral 
function  may  be  obtained  like  the  “limit”  of  a  P.r.f.  as  m  — >  +  » . 

On  the  other  hand,  from  the  above  results,  it  is  easy  to  prove  the  classical  theo¬ 
rem  of  Khintchine  which  states  a  necessary  and  sufficient  condition  that  a  func¬ 
tion  r(h )  be  a  characteristic  function  is  that  it  be  the  limit,  uniformly  over  every 
finite  interval,  of  a  function  of  the  form 

f  R(u)R(u-\-  h)  dh  , 
j  —  00 


where 
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f+C°\R(u)  \Hu=  1 
•/  —  00 


(integrals  in  the  sense  of  Lebesgue). 

2.3.  Another  application.  The  convergence  to  a  Laplace  process  may  be  easily  ob¬ 
tained  for  random  functions  of  the  following  form: 


X'{t)  =  f+C°  Y  (/,  r )  dN  (  t)  , 
—  00 


where  F(/,  r)  is  a  random  variable  (compare  1.4).  Let  us  take,  for  example,  the 
case  in  which  X'{t)  is  the  number  of  conversations  in  a  telephone  center  without 
loss,  operating  since  t  =  —  °°  (compare  1.4). 

Putting  X'{t)  =  m\  +  y/m\Z(t),  we  find  that  the  law  of  Z(t)  tends,  as 
w  — »  +  00 ,  to  a  stationary  Laplace  process,  the  correlation  function  of  which  is 

(15)  r{h)  =i/*+“[l -F(u  +  h)]du  =  l  f+C°  [1  -F(u)]du,  h>  0 

A*'  o  A  Jh 

[if  F{u)  =  1  —  £-u,  then  r(/r)  =  e~^']. 

For  h  >  0  this  function  r(/r)  is  monotonic,  nonincreasing,  lim  r(h)  =  0 

A — ►“1-00 

[ r{h )  does  nfct  vanish],  -jy  r{h)  exists. 

ah 

T'{h)  =Jh  r(k)  =  h>° 


r'(  +  0)  =  -£[l-F(  +  0)], 

r'(-h)  =  -r'(A), 


h  >  0 


For  h  >  0,  r'(^)  is  monotonic,  negative  and  r'(+  °° )  =  0.  From  these  results,  it 
is  easy  to  give  a  new  proof  of  a  theorem  of  Polya  [6,  p.  115]. 


3.  Transforms  of  a  random  function 

Suppose  that  the  random  function  X  (t)  represents  a  voltage  applied  to  the  input 
of  a  quadratic  detector.  It  is  known  that  the  response  voltage  is  then  [X(/)]2; 
therefore,  we  are  led  to  the  following  problem : 

Let  F(a:)  be  an  ordinary  function  of  x  defined  in  —  «>  <a;<+0°,  and  let  X(t) 
be  a  random  function,  not  necessarily  a  P.r.f.,  Y(t)  =  F[X(/)]  is  a  random  func¬ 
tion.  The  problem  is  to  study  Y(t),  knowing  X(t).  The  main  problem  is  generally 
to  find  the  temporal  law  of  Y{t)  from  the  one  of  X(t).  This  last  problem  has  been 
considered  by  Blanc  Lapierre  who  proposes  to  evaluate  the  moments  of  the  ran¬ 
dom  variable  F ( t )  and,  more  generally,  those  of  the  ^-dimensional  random  variable 
[F(/i)  .  .  .  F(/n)].  About  this  computation  we  make  the  following  remark:  Let  us 
suppose  that  there  exists  a  function  U(v),  the  Fourier  transform  of  F  (*),  such  that 

/-  4*oo 

F(jc)  =  /  eivxU(v)dv. 
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Let  <j>(h ,  vi\  h,  Vi\  .  .  .  ;  tnvn)  be  the  characteristic  function  of  the  ^-dimensional 
random  variable  [X(/i),  .  .  .  ,  X (/„)].  We  may  write  formally 

(16)  E[Y(h)  Y  (h)  •  •  •  Y  (4)  ] 

—E  {  V  [X  (k)  ]  X  V  [X  (/2)  ]...XF[X  (/„)  ]  } 

-fij  ...  X  [f+Je^U(vn)dvn]  | 

=  E  j  J.  )h.  f  exp  i  [  ViX  (h)  +  .  .  .  +  vnX  (tn)  ]  U(  vx)  X  .  .  .  U(  vn)  d  vx  .  .  .  d  vn  | 


=/...>«„ 


®1  ®n)  U  (  V\)  .  .  .  U  (vn)  dv!  .  .  .  dvn. 


It  is  easy  to  extend  this  result  to  the  case  in  which  the  //s  are  not  all  different. 
Then  we  have  the  formula  to  compute  any  moment  attached  to  F(/),  from  U  (v) 
and  the  characteristic  1,  vn)  of  X(t).  It  is  easy  to  give  sufficient 

conditions  to  validate  the  above  formal  computation.  For  example,  it  is  sufficient 

/+OO 

|  U(v)  |  dv  <  +  <» .  In  particular  it  is  easily  seen  that 

-CO 

E  [  Yn(t)  ]  =  f.  .  .  vu  .  .  .  ,  vn)  U(v  1) .  .  .  U  (  vn)  d  Vi  .  .  .  d  vn . 

^  n  ^ 


We  may  compare  (16)  and  the  direct  formula  (17), 

(17)  E[Y{h)..  .  Y  (4) )  ...V{xn)dF{tu  xx tn,xn), 


where  F(h,  X\  \ . . .  ;tn,  xn)  stands  for  the  distribution  function  of  the  ^-dimensional 
random  variable  .  .  .  ,  X(/n)}.  (16)  is  better  than  (17)  whenever  </>(/i,  Vi, 

.  .  .  ;  tn,  vn)  is  better  known  or  simpler  than  F(h,  Xi\  .  .  .  ;  xn),  which,  in  par¬ 
ticular,  is  the  case  if  X(t)  is  a  P.r.f.,  since  v\,  .  .  .  ;  t„,  vn)  is  then  given  by  (9). 

It  follows  that  one  may  compute  the  moments  of  some  random  variables  con¬ 
nected  with  Y(t),  particularly  those  of  the  linear  functionals  of  F(/);  for  example, 
one  has  often  to  consider  random  variables  such  as 

1  CT 

Zr  —  Tf,!  Y  (/)  dt  (stochastic  integral) , 

1  J  o’ 

the  temporal  mean  of  Y(t)  in  (0,  T ).  If  this  stochastic  integral  has  a  meaning,  it  is 
easy  to  construct  a  method  to  prove  the  result  obtained  by  Blanc  Lapierre  that 

if  R(t,  t )  depends  only  on  t  —  t,  X(t)  =  f  R(t  —  T)dN(r),  and  V(x)  satisfies 

some  wide  suitable  conditions,  then  ZT  is  asymptotically,  as  T  — »  +  °o ,  a  Laplace 
randpm  variable. 

According  to  results  of  section  2  we  conceive  that  the  preceding  results  can  in 
certain  cases  permit  the  study  of  a  random  variable  of  the  form 

(18)  H=  fhV[XAt))dt, 

where  X\ (/)  is  a  Laplace  process  of  any  kind. 
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Remark.  This  leads  us  to  make,  incidentally,  the  following  remark.  Let  V(x )  = 
x2.  Kac  and  Siegert  [4]  have  shown  that  if  X\ (t)  is  a  Laplace  stationary  process, 
the  correlation  function  of  which  is  r(h),  the  characteristic  function  <f>(v)  of  H  is 
equal  to  [D(2iv)]~l/2,  where  D(\)  is  the  Fredholm  determinant  of  the  integral 
equation  with  kernel  r(t  —  r).  But  this  result  is  not  limited  to  Laplace  stationary 
processes. 

Let  us  suppose  that  X\ (t)  be  Laplacian  and  real,  for  instance  with  covariance 
r(/,  r),  continuous  on  the  finite  interval  [a,  &].  Let  (Xi)  be  the  eigenvalues  of  the 
integral  equation 

(19)  X  (/,  T)j{T)dr  =  f{t)  +|(0, 

Ja 


and  let  [/,(/)]  be  the  normed  corresponding  eigenfunctions 


Let  us  set 


f  fi(t)  fj(t)  dt=  da . 

J  a 

Xi  =  J  fi(t)Xi(t)dt  in  quadratic  mean  . 


We  know  [7,  pp.  327-328]  that 

E[XiXi]  =0  if  i*j 

_  J_ 

MO  Mr) 


(20) 


r«,  t)  = 


X, 


if  i  =  j . 
for  a  St,  t^b] 


Xi (/)  =  ^Mt)Xi  in  quadratic  mean  for  a  S  t  S  b,  so  that 

i 

(21) 

i 

In  an  heuristic  manner  we  can  reason  in  the  following  way:  X ,  has  for  its  prob¬ 
ability  density 

V2xXi  e 

and  the  Xi  s  are  independent  of  each  other,  so  that 


*(»)=£[  e<°"i  =  f:j:  n  exp  I  ■ i  ^ **5  • -  w,  1  •  ■  •  ix‘  ■ 


the  formal  computation  of  which  is  immediate  and  gives 


(22) 


<f>  (  v) 


_ 1 _ 

VD(2iv)  ’ 


if  D(\)  is  the  Fredholm  determinant  of  the  equation  (19).  This  demonstration  can 
be  made  more  rigorous,  but  not  in  an  immediate  manner.  Another  method,  not 
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very  elegant,  but  immediately  rigorous,  is  the  following.  From  (21),  we  easily  get 
expressions,  in  terms  of  the  Xfs,  of  moments  and  then  of  cumulants  of  H;  whence, 
a  development  in  series  of  log  <f>(v )  with  the  following  form 


(23) 


log  <(>  (  v)  =  T.  vn . 


where  the  vn’ s  are  the  cumulants  of  H.  On  the  other  hand,  let  us  call  Tn(/,  r)  the 
»-th  iterated  of  the  kernel  T(/,  r)  [on  the  interval  (a,  &)].  We  know  that 


r„(/,  r)  =  £  — 


X*? 


At  last,  let  us  set 


An  =  jfVn (/,  t)dt= 


We  know  that  (cf.  p.  237  in  [8] ) 


-jr  [log  D  (\)  ]  =  -^A„X»-', 


dX 

whence  a  development  of  log  [D(2iv)~1/2\ 

(24)  log  [D  (2iv)  ~1/2]  =  ]^2 Sn~-  v* . 


The  identity  of  the  developments  (23)  and  (24)  is  easily  verified;  this  proves  the 
equality  (22).  We  have  supposed  a  and  b  finite,  but  it  is  easy  to  give  sufficient  con¬ 
ditions  in  order  that  the  result  be  yet  available  for  a  or  b  or  both  a  and  b  infinite. 
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